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ABSTRACT 



The anomalous magnetic moments of both positive and negative muon have been 
meausred to a precision of 0.7 parts per million (ppm) at the Brookhaven Alternat- 
ing Gradient Synchrotron. The results are over an order of magnitude more precise 
than the previous measurement of the muon. The results = 11 659 204(7) (5) x 
10~ 10 (0.7ppm) and a„- = 11 659 214(8)(3) x 10~ 10 (0.7 ppm), where the first un- 
certainty is statistical and the second is systematic, are consistent with each other. 
The average for the muon anomaly is a M (exp) = 11 659 208(6) x 10~ 10 (0.5 ppm). 
The standard model value of is calculated to a precision of 0.6ppm, which is dom- 
inated by uncertainty of hadronic contributions. There is a significant disagreement 
between the hadronic contribution calculated from the hadronic cross section of e + e~ 
collision and that from hadronic r decay. The difference between experimental and 
theoretical values of a M is 2.4<r for e + e~ data or 0.9cr for r decay data. 
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1 Introduction 



The gyromagnetic factor of a particle is defined as 

»=fci&s)/(f)' (1) 

where /x, S, m, and e are the magnetic moment, the spin, the mass, and the charge 
of that particle. For a pointlike lepton Dirac theory gives g = 2. Virtual radiative 
corretions cause the g value for leptons to differ from 2. The anomalous g- value is 
defined as a = (g — 2)/2. 

The anomalous (/-value of the electron and muon, a e and a M , played an 
important role in the development of particle physics and continue to serve as fun- 
damental quantities for testing the validity of the Standard Model and putting 
stringent constraints on speculative theories beyond the Standard Model. The value 
of a e is measured to the precision of 4 parts per billion (ppb). This value is currently 
the most precise measurement of fine structure constant a. 

For the muon, because the muon mass m M is larger than the electron mass 
m e , the contribution to a M from heavier particles such as hadrons or weak vector 
bosons through virtual processes is greater than that to a e typically by a factor of 
(m M /m e ) 2 ~ 4 x 10 4 . This sensitivity enhancement of a M applies generally also for 
postulated new particles such as supersymmetric particles. Hence, the measurement 
of a M provides us more information than that of a e though a M is less well known than 
a e . 

The muon anomalous magnetic moment was measured three times in CERN 
[H I210|) an d the precision of the last CERN measurement is 7.2 ppm. The muon 
g — 2 experiment at Brookhaven National Laboratory (E821) measures the a M of 
both the positive and the negative muon to 0.7 ppm. Two values agree with each 
other and the precision of the combined value of a M is 0.5 ppmjU 13 IE]- 

2 Principle and setup of the experiment 

For a muon moving in a uniform magnetic field B, which is perpendicular to the 
muon spin direction and to the plane of the orbit, and with an electric field E, both 
the spin and the momentum precess. The angular frequency difference, uj a , between 
the spin precession frequency u s and the cyclotron frequency u c is given by 

-> e 



a^B 



7 2 - 1 



(3 x E 



(2) 
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for /3 ■ B = and (3 ■ E = 0, where and e are the mass and charge of the muon, c 
is the speed of light, f3 is the ratio of speed of the muon and the speed of light, and 
7 = 1/a/1 — /3 2 . In muon g — 2 experiment, the dependence of u a on the electric 
field is eliminated by storing muons with the "magic" 7 = 29.3 such that 

a„ - = 0. (3) 

This corresponds to a muon momentum p = 3.09 GeV/c. Hence, the measurement 
of oj a and B determines a^, 

«, = (4) 

In our experiment, the magnetic field is measured by an NMR system which is 
calibrated with respect to the free proton NMR angular frequency uj p as 

B = 5Sfc, (5) 
2/x p 

where /x p is the magnetic moment of the free proton. By combination with El and 
using muon spin S/% = 1/2 and the definition of — 2(1 + o M ), 



~^c~2^ hg^2^ (l+a M )Mp /W^P 



so that 



Q M ~ P hu)„ ~ 4/x M na; p ~ Mm _ 7 /.. V 1 i" °W> [[) > 

(l+a M ) 

<V = / ]—■ ( 7 ) 

The ratio fi^/fi p = 3.183 345 39(10) [Qj, as determined from muonium spectrocopy 

ma. 

The frequency uj a is measured by counting high energy decay positrons. In 
the muon rest frame, parity violation in the decay [i + — > e + u e i? tl causes e + to be 
emitted preferentially along the muon spin direction. Positrons emitted along the 
muon momentum direction get the largest Lorentz boost and have highest energy 
in the laboratory frame. Hence, the angle between the muon momentum and muon 
spin determines the number of high energy decay positrons. While the muon spin 
precesses with angular frequency uj s and the muon momentum precesses with angular 
uj c in the magnetic field, the number of decay positrons with energy greater than 
some threshold energy E oscillates with the angular frequency u a = oj s — uj c and is 
ideally described by 

N(t) = N (E)e~ t ^{l + A(E) sm[u a t + ME)}} (8) 
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in which 7T is the dilated muon life time, A is the asymmetry, <fi is the phase of the 
oscillation. Both A and <fi depend on the energy threshold of the decay positrons. 
However, the frequency uj a does not depend on E. 

The general arrangement of the muon g — 2 experiment is shown in Fig. ^ 
The proton beam from the Alternating Gradient Synchrotron(AGS) strikes a nickel 
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Figure 1: General arrangement of the muon g — 2 experiment. 

target to produce pions. The pions are captured into a secondary beamline where 
a fraction of about 50% decays to muons. We select the high energy muons which 
come from the forward decay. They are polarized due to the parity violation in the 
decay n + — > /i + z/ M . The muons are injected through a superconducting inflector [TT] 
into the muon storage ring with 1.45 T uniform magnetic field [T2] which is measured 
by a NMR system [TB*] relative to the free proton NMR frequency. A pulsed magnetic 
kicker Jl] at about 90° from the injection point gives muons a 10 mrad deflection 
which places the muons onto stored orbits. The electrostatic quadrupole[T5] fields 
provide vertical focusing. The decay positrons are detected by 24 lead/scintillating 
fiber electromagnetic caloriment ers 1 1 6] read out by waveform digitizers. Both the 
waveform digitizer and NMR clocks were phase locked to the Loran C frequency 
signal. 



3 Data Analysis 

The data analysis in the experiment is divided into two independent parts, u p and u a 
analyses. The absolute values of ui p and u) a from analyses were modified by arbitrary 



5 



offsets. No one knew both absolute values of u p and u a . The values of R = u a /u p 
and then a M were evalueated only after both analyses had been finalized. 

3.1 u p analysis 

The measurements of the magnetic field at sub-ppm lever is done with NMR probes 
and is expressed relative to the free proton NMR frequency u p . 

The NMR frequencies are measured with respect to a standard probe, 
which is a spherical water sample. The calibration of the standard probe with 
respect to the free proton is obtained from other experiments fTT] . 

A trolley with 17 NMR probes was used to measure the field around the 
storage ring. The 17 trolley probes are calibrated with respect to the standard 
probe through a secondary probe called plunging probe, which can be mounted 
on a stand and moved in radial and vertical direction in vacuum chamber at one 
azimuthal location. The calibrations were done by comparing the NMR readings of 
the plunging probe and the trolley probes taken at the same location. The errors 
come from the position uncertainties of the trollye probes and of the plunging probe, 
and from the B field inhomogeneity, Large magnetic sextupole moments and field 
gradient in azimuth direction were introduced to measure the relative position of 
active volume of trolley probes and the plunging probe. The effect on the trolley 
probe calibration from temperature, power supply and magnetic field were also 
carefully studied. 

We measured the magnetic field with trolley about twice per week. One 
of the measurements of the magnetic field with the trolley center probe is shown vs. 
azimuth in Figure 121 For most of the ring, the magnetic field varies about ±50 ppm. 
The field measurements are averaged over azimuth for all 17 probes individually. 
Since the field component in radial and azimuthal direction are much smaller than 
that in vertical direction, the field can be expanded in 2-dimension 

B(x, y) = B(r, 9) = B + £ a, (-Y cos(i0) + £ h (-Y sm(jd), (9) 

i=i V? V j= i \r / 

where x is the radial direction in our ring, and (r, 9) are polar coordinates. The 
center of the storage region is at r = and the direction pointing outward is 9 = 0. 
The value of r is set to be 4.5 cm. A half ppm contour plot and the multipoles of 
the azimuthal average of the field measured in one trolley run are shown in Figure El 

The field between trolley runs was measured by a fixed probe system of 378 
NMR probes. More than 120 probes are used in the analysis because of their good 
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Multipoles (ppm) 
normal skew 
Quad -0.29 0.15 
Sext -0.71 -0.48 
Octu 0.06 0.01 
Decu 1.08 0.39 



6=180° 

Figure 2: The left plot shows the magnetic field measured with the trolley center 
probe. And the right plot shows the half ppm contour plot and the multipoles of the 
azimuthal average of the field measured in a trolley run. 

qualities. The average measurements from those fixed probes during the trolley run 
is calibrated with respect to the trolley measurements. 

The muon distribution is taken into account by using the field value at 
the center of the muon distribution. The radial muon distribution is measured by 
fast rotation with uncertainty of 1 mm. The vertical muon distribution is symetric 
around the central plane determined by electric quadruples, and the uncertainty is 
2 mm. 

The contributions to the total systematic error from different sources are 
given in Tabled for 2001 run. And the result for u p weighted by the muon distri- 
bution in 2001 is found to be 

u p /(2tt) = 61 791 595(15) Hz (0.17ppm) (10) 
Table 1: Systematic uncertainties for u p analysis of the 2001 data 



Source of errors 


Size [ppm] 


Absolute calibration of standard probe 


0.05 


Calibration of trolley probes 


0.09 


Trolley measurements of Bq 


0.05 


Interpolation with fixed probes 


0.07 


Uncertainty from muon distribution 


0.03 


Others ' 


0.10 


Total Systematic Error on uj p 


0.17 



Higher multipoles, trolley temperature and voltage response, and eddy currents 
from the kickers, and time- varying stray fields.. 
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3.2 oo a analysis 



In 2001 run, We collected about 4 billions of electrons 32/zs after the injection and 
with reconstructed energy larger than 1.8GeV. The detection time of electrons are 
randomized by cyclotron period to eliminate the bunch structure at early time. And 
the time spectrum after randomization is shown in Fig. El 
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Figure 3: The time spectrum of electrons 

The leading characteristics of the time spectrum are those of muon decay 
and the muon g-2 oscillation, u a , as shown in Eq. ^2 

N(t) = N e- thT [l + Acos{uj a t + 0)] (11) 

Additional effects includes pulse overlap, detector gain stability, muon loss and co- 
herent betatron oscillation (CBO). Among those, CBO in horizontal plane has the 
largest effect on the determination of u> a . CBO is caused by injecting muons through 
narrow aperture of the inflector into larger aperture of the strorage ring. Since the 
calorimeter acceptance varies with the decay position of muons, the time and energy 
spectra of decay electrons are modulated with the CBO frequency. The beating fre- 
quency between CBO and u a is very close to uj a . Hence, the CBO has large effect on 
determination of u a . Since the phase of CBO oscillation varies from to 2tt around 
the ring, its effect in the sum of spectra from all detectors is significantly smaller 
than that in any individual detector. The CBO effect is also accounted in the fitting 
functions. In 2001, we ran our experiment at two values of CBO frequency, one with 
CBO frequency larger than 2 x u a and the other with CBO frequecy smaller than 
2 x uj a . Fitting the sum spectrum with function taking into account the CBO gave 
small systematic error for both cases. 
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Two approaches were used to obtain the value of u a . In one approach, we 
fit the time spectrum directly with some multi-parameter function. The function 
includes Eq. ^2 and the effects from gain change, CBO, muon losses. The fitting 
function is slightly different for different analyzer. The other approach is called 
ratio method. The data are randomly assigned to four statistically independent 
subsets rii to 714. The subsets are rejoined in u(t) = n\{t) + n 2 (t) and v(t) = 
n 3{t ~ r /2) + n±(t + t/2), where r is an estimate of the (g — 2) period, and then 
combined to form the time spectrum r(t) = [u(t) — v(t)]/[u(t) + v(t)}. The {g — 2) 
rate modulation of v is 180° degrees out of phase compared to that of u, and to 
sufficient precision r(t) can be described by 

r(t) = A sm(uj a t + <j) a ) + (r a /16r) 2 . (12) 

where r a is the (g — 2) period and r is the dilated lifetime of muons. The ratio 
r(t) is largely insensitive to change of observed counts on time scales larger than r . 
However, CBO effect still has to be taken into account in ratio method. 

We completed four independent analyses. They all agree with one another 
within the expected variation. The combined result 

uj a = 229 073.59(15)(5)Hz (0.7ppm), (13) 

where the first uncertainty is statistical uncertainty and the second is the combined 
systematic uncertainty. And the systematic errors are listed in Table |21 

Table 2: Systematic uncertainties for uj a analysis of the 2001 data 



Source of errors 


Size [ppm] 


Coherent betatron oscillation 


0.07 


Pileup 


0.08 


Gain change 


0.12 


Lost muons 


0.09 


Others t 


0.11 


Total Systematic Error on uo a 


0.21 



t AGS background, timing shifts, E field and vertical osci 



ations, beam 



debunching/randomization, binning and fitting procedure. 



3.3 results 

The value of a M was determined after the analyses of uj p and u a had been finalized, 



a„- = 11 659 214(8)(3) x 10~ 1U (0.7 ppm) 



(14) 
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in which the correlations between the systematic uncertainties have been taken into 
account. 

In Eq. fM^/fip is measured by [TUj to high precision for fi + only. Hence 
a M + and a M - cannot be done directly. Instead, we compared R = uj a /uj p values from 
fi + ad n~. The difference is 



e = - 



2 + Rp- 



-9.4 ±9.7 x 10" 



(15) 



and is consistent with 0. Assuming CPT, we combine the results from /i + and fi 
and obtain 

a M = 11 659 208(6) x 10~ 10 (16) 
Fig. E] shows the measurements and the theoretical prediction of a M . 
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Figure 4: The experimental results and the theoretical prediction of a 



i-i- 



4 Theory 

In the Standard Model, the contributions to a p are divided into three types, QED, 
hadronic and weak. 



a^SM) = a^(QED) + a^(had) + a M (weak). 



(17) 



The QED contribution has the largest contribution but least uncertainty. 
The calculation up to 5 th order of a/(2n) gives 



a M (QED) = 11 658 470.56(29) x 10 



-10 
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A recent recalculation of the terms of 4 th order of a/(2ir) increases the QED contri- 
bution to 

a M (QED) = 11 658 471.94(14) x Kr 10 [E] (19) 
The weak contribution has been calculated to the second order and is evaluated as 



a M (weak) = 15.4(2) x 10 1U [2 



(20) 



The hadronic contribution cannot be calculated from pertubative QCD 
alone because it involves low energy scales near the muon mass. However, the 
leading order of hadronic contribution a M (had, lo) (2TJ, |22l EHj can be determined 
directly from the annihilation of e + e~ to hadrons through a dispersion integal 



a,(hadlo) = (^) 2 f° ^K(s)R(s 

V 3ll ) JAml s 2 



(21) 



in which s is the energy in the center of mass frame, K(s) is a kinematic factor and 



R(s) 



^totaiie^e 



hadrons) 



l 1 I 1 ' 



(22) 



The value R{s) can be indirectly determined using data from hadronic r decays, 
the conserved vector current hypothesis, plus the appropriate isospin correction. 
In principle, the large statistics of r decay data could improve the precision of 
a M (had). However, discrepancies between the r and e + e~ data exist as shown in 
Fig. El The e + e~ data were dominated by the data from CMD2 experiment. The 
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Figure 5: The discrepancies between r and e + e~ data at different energy scales. The 
data points are e + e~ data from different experiments, and the green bands are the r 
data. 



recent result from radiative return data[22] with comparable uncertainty agree with 

11 



the CMD2 experiment. The next leading order hadronic contribution [26j 1231221 and 
the hadronic light-by-light scattering [23 I2H1 1201 HH 122] are calculated by serveral 
groups and the combined the results are 

a^had nlo) = 10.0(0.6) x 10~ 10 a^lbl) = 12.0(3.5) x 10~ 10 (23) 

The large discrepancies between r and e + e~ data prevent us from com- 
bining two to obtain a single standard model prediction. The differences between 
experimental measurement and the standard model prediction are 

Aa^e+e") = 23.9(7.2 had i )(3.5 lbl )(6 cxp ) x lO" 10 (2.4a) (24) 

Aa,(r) = 7.6(5.8 had i )(3.5 lbl )(6 exp ) x lO" 10 (0.9a) (25) 

The question concern the discrepancies betweent e + e~ and r data are under further 
theoretical scrutiny and we expected more radiative return measurements and efforts 
from latice calculation. 
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